Origin of low-temperature magnetic ordering in Gai-a^Mn^^N 
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By employing highly sensitive millikelvin SQUID magnetometry, the magnitude of the Curie 
temperature as a function of the Mn concentration x is determined for thoroughly characterized 
Gai-ccMuccN. The interpretation of the results in the frame of tight binding theory and of Monte 
Carlo simulations, allows us to assign the spin interaction to ferromagnetic superexchange and to 
benchmark the accuracy of state-of-the-art ah initio methods in predicting the magnetic character- 
istics of dilute magnetic insulators. 



The extensive studies of dilute magnetic semiconduc- 
tors (DMSs) and oxides over the last decade [1, 2] have 
persistently confronted the researchers with experimen- 
tal and conceptual challenges, making the field to be one 
of the most controversial in today's condensed matter 
physics. More specifically, it has become increasingly 
clear that the premise of dilute magnetic alloys - where 
the magnetic constituents incorporate randomly and sub- 
stitutionally into the host crystal - breaks entirely down 
in a number of systems [3]. In particular, the distribu- 
tion of magnetic ions is often non-uniform and the ions 
tend to occupy also interstitial positions. Moreover, the 
paramount importance of disorder, defects, and strong 
correlation makes that the theoretical and computational 
modeling of these materials has often been misleading. 
However, recently some consensus on the ah initio tools 
appropriate to study the (ferro) magnetism in these sys- 
tems has been reached [1, Q . 

Since GaN and its alloys with Al and In have al- 
ready realized their potential in photonics and high 
power electronics, reaching the status of the technolog- 
ically most significant semiconductor materials next to 
Si, the addition of magnetism opens wide application 
prospects. In particular, the presence of ferromagnetic 
interactions withouthd^nd carriers, together with a sizable 
spin splitting of the excitonic states revealed already for 
Gai-a^Mn^^N 0,0] 7 indicates the suitability of this sys- 
tem for magnetooptical devices, such as optical isolators, 
circumventing the destructive effect of antiferromagnetic 
interactions specific to II- VI Mn-based DMSs [sj. Sur- 
prisingly, however, in previous works a variety of different 
magnetic behaviors is reported for Gai-^^Mn^^N at the 
same nominal Mn concentration x. This compound was 
found by some groups to be non-magnetic [9], whereas 
according to others it shows either low-temperature spin- 
glass freezing [13] or ferromametism with a Curie tem- 
perature Tc ranging from 8 K jllj up to over 300 K [2l.[l2|. 



In this Letter we report on studies of magnetic hystere- 
sis down to millikelvin temperatures for Gai-a^Mua^N epi- 
taxial layers, in which the high crystallinity, the random 
distribution of Mn ions, and the extremely weak degree of 
compensation by residual donors were assessed by a range 
of electron microscopy, synchrotron radiation, optical, 
and magnetic resonance techniques [6]. The magnetic 
phase diagram Tq {x) established in this way allows us to 
verify the predictive power of ah initio methods. With 
the support of tight-binding theory and Monte Carlo sim- 
ulations, we corroborate the experimental results and 
find that state-of-the-art first principles approaches over- 
estimate the magnitude of the Mn-Mn exchange energies 
by an order of magnitude. In this way, we signify that 
dilute magnetic insulators constitute a relevant system 
to benchmark newly developed tools for computational 
design of functional magnetic materials. 

The samples discussed here have been grown by met- 
alorganic vapor phase epitaxy according to the procedure 
reported previously 0,0, [l^. In particular, Gai-a^Mn^^N 
has been deposited onto GaN/c-sapphire at a substrate 
temperature of 850° C. In order to maximize the homo- 
geneous and substitutional incorporation of Mn, to avoid 
phase separation, and at the same time to vary from sam- 
ple to sample the actual concentration of Mn in a con- 
trolled way, the flow rate of the Ga precursor (TMGa) 
has been changed over the samples series from 5 to 1 
standard cubic centimeters per minute (seem), the tem- 
perature of the Mn precursor source (MeCp2Mn) from 
17°C to 22° C, while its flow rate has been maintained 
constant at 490 seem for all the samples considered. 

The films have been thoroughly characterized by 
secondary-ion mass spectroscopy; high resolution (scan- 
ning) transmission electron microscopy with capabili- 
ties allowing for chemical analysis, including energy- 
dispersive x-ray spectroscopy, high angle annular dark- 
field mode, and electron energy loss spectroscopy; high- 
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FIG. 1. (Color online) (a) Sketch of the mK magnetometer employed for this work. For maximum linearity of the set-up 
the superconducting pick-up coils are made of pure niobium wire and the magnetic field is generated by a copper wound 
electromagnet immersed in a helium dewar. An additional copper-made flux-balancing coil (c. c.) is placed asymmetrically 
with respect to the main coil and is connected to the same power supply through a resistive current divider (Rl and R2) adjusted 
to obtain an adequately flat magnetic field (reference) response when a reference undoped GaN layer is placed inside the pick- 
up. The whole set-up has been calibrated against the known magnetic moment of a signal coil located at the sample position. 
(b)-(d) Magnetic hysteresis loops at various temperatures for Gai-ccMn^N measured in the dilution fridge set-up (lines) and in 
a commercial SQUID magnetometer (bullets), respectively. The reported values of the magnetization are normalized to their 
magnitudes at 70kOe and at 1.85 K. 



resolution and synchrotron x-ray diffraction; synchrotron 
extended x-ray absorption fine-structure; synchrotron x- 
ray absorption near-edge structure; infrared optics; elec- 
tron spin resonance; and superconducting quantum inter- 
ference device (SQUID) [6]. This extensive analysis has 
allowed us to rule out the presence of Mn-precipitation 
and indicates that: (i) up to a Mn content of 3.1% at 
least 95% of the Mn ions has the charge state 3+ and 
is substitutionally incorporated in the host crystal; (ii) 
the layers are highly resistive even at 300 K indicating 
no charge transport via band or Mn gap states. The con- 
centration X of Mn has been tuned from 0.5% of cations 
(2x10^^ cm~^) to 3.1% for different samples. Here, the 
layers with respectively x = l.l% (sample 1080), 1.8% 
(sample 1106) and 3.1% (sample 1142) are considered, 
where the code of the sample numbers is the same 
adopted previously Q. 

Previous magnetization studies of these Gai-a^Mn^^N 
films were carried out down to 1.8 K [6| , employing a com- 
mercial SQUID magnetometer and a measurement tech- 
nique specifically developed in order to examine mean- 
ingfully thin layers of magnetically dilute semiconductors 
14j . The data pointed to the presence of ferromagnetic 
interactions between Mn spins but no signatures of long 
range magnetic order were found In order to extend 
the range of our measurements, we have installed in our 
^He/^He dilution fridge a homemade SQUID-based mag- 
netometry set-up that allows for high-sensitivity detec- 
tion of magnetic moments during continuous field sweeps 
up to about 800 Oe and down to 20 mK, as sketched in 



Fig.IHa). 

In Fig. [Hb)-(d) the magnetization loops obtained at 
various temperatures for the studied films are reported. 
The emergence of a long-range ferromagnetic order is wit- 
nessed by the appearance of hystereses, whose width (and 
height) exhibit a critical behavior on lowering tempera- 
ture, allowing to determine the magnitude of Tq for par- 
ticular films, as shown in Fig.O In contrast to the case of 
those dilute magnetic semiconductors and oxides which 
show high values of Tq independent of the magnetic ion 
concentration x [2], our data point to a significant vari- 
ation of Tq with X, as expected from any model of fer- 
romagnetism within a system of randomly distributed 
localized spins 0, [HI . 

To put these data into a broader context, we note that 
over the last decade Gai-a^Mn^^N has reached the status 
of a model system, whose magnetic properties have been 
theoretically examined by over twenty different groups 
employing a variety of ab initio methods [1, [H, [l6|, [l7j . 
Notably, these approaches predict consistently the pres- 
ence of ferromagnetic coupling, as observed. Moreover, 
again in agreement with experimental observations, a 
cross-over to antiferromagnetic interactions is expected 
upon donor compensation that leads to the reduction 
Mn^+ Mn^+ |16j. However, a quantitative comparison 
of the experimental and theoretical Tq values as shown 
in Fig. [2] implies that the state-of-the-art computational 
approaches overestimate the measured values by an order 
of magnitude. 

To unravel this issue, we recall that uncompensated 
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FIG. 2. (Color online) Magnitude of the coercive field as a 
function of temperature for the considered Gai-ccMn^N sam- 
ples (full points), employed to determine the Curie tempera- 
tures. Inset: experimental Curie temperatures as a function of 
Mn content x, together with the experimental result of Ref.fTll 
(solid squares and circle, respectively). The dotted line shows 
the scaling dependence Tc oc with m = 2.2. The results 
of ab initio (Ref. 4, open diamonds) and tight-binding (this 
work, open squares) approaches are also shown. 



Gai_a;Mna;N can be classified as a dilute magnetic in- 
sulator 0, in which the absence of electrons and holes 
makes carrier-mediated spin-spin coupling [18] irrelevant, 
and the lack of mixed valence - all magnetic ions are in 
the same 3+ charge state - precludes the presence of 
double exchange [l9[. In this situation, superexchange 
accounts for spin-spin interactions [20|. Its sign is de- 
termined by the Anderson-Goodenough-Kanamori rules, 
whereas the character of the p — d hybridization controls 
the magnitude of the scaHng exponent m describing the 
dependence of the magnetic critical temperature on x, 
Tc{x) (X x^. The fact that the scaling law with similar 
exponents is obeyed by spin glass freezing in II- VI DMSs 
(m = 1.9 ± 0.1) [21] and by ferromagnetic ordering in 
Gai_^Mn^N (m = 2.2 ± 0.3) - as evidenced in Fig. [2]- 
strongly supports the superexchange scenario. 

In order to evaluate the sign and magnitude of the spin- 
spin coupling, we adopt for Gai-a^Mua^N an experimen- 
tally constrained procedure developed by one of us and 
co-workers for II- VI compounds doped with transition 
metals (TM) [22[. Within this approach, the magnetic 
ions are described in terms of the Parmenter's [23| gen- 
eralization of the Anderson hamiltonian for the relevant 
electronic configuration of the TM taking into account 
the Jahn- Teller distortion 24, 25|, whereas the host band 
structure is modeled by the tight binding approx- 

imation, employing the established parametrization for 
GaN [26| in the cubic approximation. Other parameters 
of the model [22[ are taken from experimental studies of 



optical [27| as well as photoemission and soft x-ray ab- 
sorption spectroscopy [28| of Gai-^^Mn^^N. In particular, 
the charge transfer energy between the Mn ion and the 
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FIG. 3. (Color online) Exchange energies provided by the 
tight binding model for Mn pairs in zinc blende GaN as a 
function of the distance between Mn spins in units of lattice 
parameter a. The numbers of equivalent cation sites at the 
particular distances are also given. 



top of the valence band, Mn^+ Mn^+ is ei = -1.8eV 
[27|, which together with the on site correlation energy 
for Mn^+ ions [27], U = 1 eV and the on site exchange 
energy for Mn^^ ions, A = E{S = 5/2) - E{S = 3/2) = 
2eV leads to 62 = 4.8 eV, where the uncertainty on the 
relevant energies, ei and 62, is presumably of the order of 
±0.5eV. The magnitude of the p-d hybridization energy 
is Vpda = -1.5±0.1eV [28;]. 

Since the effect of spin-orbit splitting is small in the 
valence band of GaN, the spin dependent interaction be- 
tween two Ga-substitutional Mn spins is described by 



a scalar Heisenberg coupling H- 



7(5 



where 



7 and S denote the one t2g orbital (either xy^ xz or yz) 
which is empty at the Mn^+ ion i and j , respectively. The 
magnitudes of J^-^ are evaluated within the fourth order 
perturbation theory in Vpd for all possible orbital config- 
urations 7 and 5. Similarly to the case of Cr^+ ions in 
II- VI compounds [22|, the main contribution originates 



from quantum hopping involving occupied t2g orbitals at 
the one Mn^+ ion and the empty orbital at the other 
Mn^+ ion. For the orbital configurations in question we 
find that the interaction is ferromagnetic at all distances. 

In order to compare quantitatively the theoretical and 
experimental results, we assume a statistical distribution 
of directions corresponding to tetragonal Jahn- Teller dis- 
tortions and determine an average value of the exchange 
energy Jij characterizing the coupling of Mn^+ pairs at a 
given distance i^ij in the fee cation sublattice. In this way 
we obtain the values of Jij shown in Fig. [3l By taking 
into account the coupling up to ten subsequent neigh- 
bor positions we allow for the formation of a percolation 
cluster down to x ^1.2% [29|. Iniportantly, when com- 
pared to the ab initio results 0, [isl, our values of 
Jij are significantly smaller and may, therefore, lead to 
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a better agreement with the experimental data. This is 
substantiated by the Tq values summarized in the inset 
to Fig. [21 which have been obtained by employing Monte 
Carlo simulations and the cumulant crossing method (3Q| . 
We note here, that the magnitudes of are rather sen- 
sitive to the input parameters. For instance, by changing 
62 from 4.8 eV to 4.4 eV, i. e., within its expected un- 
certainty, the computed Tq values are in agreement with 
the experimental data. 

In summary, the substantial agreement between our 
experimental values of Tc{x) and our tight binding 
and Monte Carlo simulations, has made it possible to 
identify ferromagnetic superexchange as the microscopic 
mechanism accounting for the ferromagnetic interaction 
between localized spins in Gai-a^Mn^^N. Because of 
its short range character, this coupling leads to rather 
low Tc{x) values. Furthermore, the results allow to 
shed light on the predictive power of the current first 
principles methods, broadly employed to treat the case 
of Gai-a^Mn^^N. We note, that the computed magnitudes 
of Tc(x), obtained from the available implementations 
of DFT, are much higher than the experimental val- 
ues. This disagreement originates presumably from 
an overestimation - inherent to the local spin density 
approximation - of the metallization of the d level, this 
being an error not affecting the approach employed 
here. It would be now worth to assess, whether other ab 
initio methods, like e. hybrid density functionals [25| 
or a combination of DFT with a dynamic mean field 
approximation (sij, could bring computational results 
even closer to the experimental values. 
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